Emerging non-equilibrium bound state in spin current-local spin scattering 
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Magnetization reversal is a well-studied problem with obvious applicability in computer hard- 
drives. One can accomplish a magnetization reversal in at least one of two ways: application 
of a magnetic field, or through a spin current. The latter is more amenable to a fully quantum 
mechanical analysis. We formulate and solve the problem whereby a spin current interacts with a 
ferromagnetic Heisenberg spin chain, to eventually reverse the magnetization of the chain. Spin- 
flips are accomplished through both elastic and inelastic scattering. A consequence of the inelastic 
scattering channel, when it is no longer energetically possible, is the occurrence of a new entity: 
a non- equilibrium bound state (NEBS), which is an emergent property of the coupled local plus 
itinerant spin system. For certain definite parameter values the itinerant spin lingers near the local 
spins for some time, before eventually leaking out as an outwardly diffusing state. This phenomenon 
results in novel spin-fiip dynamics and filtering properties for this type of system. 
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INTRODUCTION 

Most current computer hard drives utilize a technol- 
ogy for memory storage which requires a switching of 
states involving magnetized spin. This switching is ac- 
complished through the application of magnetic fields in 
appropriate directions. A theoretical understanding of 
this process is attained reasonably well through a classi- 
cal description via the Landau-Lifshitz-Gilbert equations 
These equations constitute a phenomenological de- 
scription, since the required darnping, whose analytical 
form is even under some debate [3| , has various possible 
origins. 

Just over a decade ago, however, theoretical propos- 
als were made to accomplish magnetization switching 
through spin transfer from applied spin currents to mag- 
netized spins A semi-classical description was used: 
the spin current was described by a plane wave, while the 
magnetized thin film that was to be flipped was described 
through a classical magnetization vector. This problem 
became known as the 'spin-torque' problem; the incom- 
ing spin current exerts a torque on the local magnetiza- 
tion. It is noteworthy that in this problem a phenomeno- 
logical damping mechanism is not required to torque the 
magnetization in the direction of the incoming spin cur- 
rent — whereas the use of a magnetic field leads only 
to precession unless some damping mechanism is intro- 
duced. The experimental observation of the 'spin-torque' 
effect has met with some limited success [1, 0, H, 0] . 

Recently, a direct measurement of the spin torque vec- 
tor depending on the voltage has been made [lO|. Fur- 
thermore, the results of this experiment imply that in- 
elastic tunneling has an important effect on the spin 
transfer torque. In fact, it appears that inelastic pro- 
cesses in the spin-flip scattering are inherent 11| for fer- 
romagnetic systems. In order to realize practical ap- 
plications of the spin torque phenomenon, it is impor- 
tant to reduce the critical current required to reverse the 



magnetization of ferromagnets. A couple of experiments 
[13, have demonstrated experimental methodologies 
to decrease the critical current. As another signature of 
spin transfer, spin-torque induced magnetic vortex phe- 
nomena are also observed 14 , 1^ 1^ . 

The semi-classi cal pict ure seems to work well in a prac- 
tical sense 
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2l|. However, especially from a 



theoretical point of view, some aspects are missing. Ulti- 
mately, spin transfer is a quantum mechanical scattering 
problem, generally inelastic, and so one would like to 
understand the spin transfer process in terms of excita- 
tions of the ferromagnet. Moreover, recent experimental 
work [22] has focussed on the impact of a spin current on 
cobalt nanoparticles with diameter less than 5nm, which 
can be used to examine the spin torque exerted on iso- 
lated nanoparticles. It has also been shown that it is 
experimentally feasible to manufacture ma gne tic nanos- 
tructures (chains of 2-10 coupled atoms) [23|. In this 
case, only a fully quantum mechanical description will 
suffice because the quantum nature of the spin operator 
representing the stationary spins in the nanoparticle is 
significant. 

The scenario of an incoming (electron) spin, often mod- 
eled as a wave packet, whose spin degree of freedom is 
coupled with local spins, has been advanced by a number 
of workers [H, [H H, [13, H [13 . The coupling between 
the incoming spin and the local spins is Kondo-like, while 
the local spins are themselves ferromagnetically coupled 
via a Heisenberg exchange interaction. The model Hamil- 
tonian is 



H 



E 

<i,j><J 



Jo(Ti-Se 



N,-l 



JiSi ■ S 



e+i 



(1) 



where cj^ creates an electron with spin a at site i, Si 
is a localized spin operator at site £, and to is the hop- 
ping amplitude between nearest neighbor sites. The first 
term allows an electron (of either spin) to propagate in 
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a band that covers all space (here in one dimension), 
while the second term is responsible for the Kondo-like 
interaction between the electron and the local spins, with 
coupling constant Jq. This takes place over a finite chain 
of length Ns- Finally, the last term models the Heisen- 
berg exchange interaction with strength Ji between the 
local spins. For a ferromagnetic chain, Ji > 0. Note, 
moreover, that if so desired, both Jq and Ji can depend 
on the position of the local spin within the finite chain. 
Fig. [T] shows a schematic of this model. 

The use of a wave packet to describe the incoming spin 
degree of freedom, and the subsequent 'real-time' analysis 
of the scattering process allows us to examine the entire 
scattering process with very fine spatial and temporal 
resolution. While the present-day experimental capabili- 
ties do not quite match this fine resolution, we anticipate 
that probing on the time and length scales we use will 
be accessible in the near future. In particular, in this 
work we identify an new feature which he have identi- 
fied as a 'non-equihbrium bound state' (NEBS), whose 
characteristics would require careful experimental identi- 
fication. This phenomenon results because of an inelastic 
scattering process that is suppressed due to energy con- 
servation. While an analytical approach does reveal some 
of the properties of a NEBS, the numerical wave packet 
calculations really allow us to see the non-equilibrium as- 
pect of this phenomenon. Both calculations are presented 
here. 

This paper is organized as follows. In the following 
section we outline means by which we solve the time- 
dependent problem. Some of our earlier work 
used straightforward expansions in the basis states span- 
ning the product Hilbert space of electrons moving on a 
lattice and stationary spins confined to a small portion 
of that same lattice. The present work uses a different 
method; the exponentiated Hamiltonian operator is ex- 
panded in a series utilizing Chebyshev polynomials 31 1 . 
This allows us to easily generate large scale numerical 
results, as described in the subsequent section. We for- 
mulate the problem for an arbitrary number of stationary 
spins (in principle, representing a magnetized thin film, 
whose magnetization is being fiipped), but focus on two 
interacting stationary spins. This allows us to focus on 
the characteristic features of the larger system, including 
the NEBS, without the considerable complexity gener- 
ated by the many scattering channels present when more 
than two stationary spins are used. Snapshots of the 
propagating wave packet reveal that in a particular re- 
gion of parameter space part of the wave packet 'lingers' 
near the stationary spins. This feature is a signature of 
the NEBS. 

In the fourth section we develop an analytical approxi- 
mation to describe the same scattering process in the con- 
tinuum limit. A preliminary decomposition of the prob- 
lem, into less familiar but more useful basis states, first 
allows us to readdress the numerical results of Section 



III. This analysis identifies the NEBS with the position- 
dependent amplitude of one of these basis states. We fur- 
ther develop the analytical approximation to derive this 
amplitude, along with expectation values for the amount 
of spin-flip expected. Thus, while we lose the trans- 
parency of the time dependent (i.e. non-equilibrium) as- 
pect of the problem, we clarify some of the physics of 
the bound state part. In Section V we conclude with 
some discussion concerning experimental observation of 
this NEBS. 



THEORY 

We adopt the most straightforward approach to the 
scattering problem, and study the time evolution of a 
wave packet, defined, at i = 0, as 

ip{x) = -^=e*'=(^-=^'')e (2) 



The calculation can take several routes at this stage. 
Consistent with the tight-binding formulation, Eq. H]), 
one can define a Hilbert space (with either open or peri- 
odic boundary conditions left of the wave packet and far 
to the right of the local spins), with typically hundred's 
of lattice sites on which the itinerant spin (hereafter re- 
ferred to as the electron, or electron spin) can hop (see 
Fig. 1 for a schematic). One can diagonalize Eq. ([T]) 
on this Hilbert space and find the complete spectrum 
of eigenstates and eigenvalues with which one can con- 
struct the time evolution of the wave packet 2^ 27, 13] ■ 
However, we find that the parameter regime and maxi- 
mum possible size of the local spin chain, for example, is 
severely restricted by computational expense within this 
approach. 

Instead we choose to solve the time dependence di- 
rectly, using the formal solution 



-iHt 



(3) 



A practical implementation of this solution is through 
the series expansion 
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where an are the coefficients of a complete orthonor- 
mal set of functions denoted by F„. A very use- 
ful basis is provided by the Chebyshev polynomials, 
Tn{x) = cos(ncos-ia;), with ro(^) = ^,Ti{X) = X 
and Tn{X) = 2XTn^i{x) ~ T„_2(X) [3l|. For this ex- 
pansion to be useful, the argument X (here, a matrix) is 
required to have norm less than unity, so a scaled version 
of the Hamiltonian is required (accompanied by a scaled 
time variable): 
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FIG. 1: (Color online) A schematic of a lattice, on which an 
itinerant spin can hop (with hopping parameter to); it can 
interact with two stationary spins (indicated by downward 
pointing (red) arrows) with coupling strength Jq. The two 
stationary spins can interact with one another, with coupling 
strength Ji. 



where y = 5t, and x 
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There are two reasons for choosing this particular ba- 
sis. First, the coefficients an{y) can be written simply as 
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where the Jn{y) are Bessel functions of the first kind. 
Second, these polynomials have a recursion relation that 
allows us to use a more compact calculation of the ex- 
pansion of the exponential of the Hamiltonian, 



Using this equation we can rewrite the expansion up to 
a given order, N"^ as 
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and the matrix elements A{i,j) are defined by 

( A{i-lJ) + 2*A{i-lJ -1) mod{i- j,2) = 
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with A(0, 0) = 1. 



This formulation allows for an efficient evaluation of the 
time evolution of the wave function, such that large lat- 
tices can be studied, both for the electron spin, and for 
the stationary spin chain. 

NUMERICAL RESULTS 
(a) non-interacting stationary spins 

The result of a typical calculation is illustrated in Fig. 
[21 Here, we have used 1600 lattice sites, and, at i = 
we have 'launched' a wave packet centred around site 700 
with a width given by a = 30. The unit of length is the 



lattice spacing, which we take to be unity for convenience. 
In all our figures we also take to = 1 as our energy scale. 
All our results will utilize an initial electron wave vec- 
tor k = 7r/2, so that no wave packet broadening occurs 
[iot . The incoming electron spin has S = 1/2, and, in 
the calculations in this paper, the stationary spins have 
S' = 1/2. A series of snapshots is shown as time pro- 
gresses forward. Initially only the incoming electron is 
present, represented as a Gaussian wave packet with only 
a spin-up component (shown as a solid (red) curve for the 
first time slice at the bottom) . The initial conditions are 
such that all stationary spins (not shown but situated at 
sites 800 and 801) have Sz = —1/2 and the incoming elec- 
tron spin has Sz = 1/2. As time advances the electron 
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FIG. 2: (Color online) Time evolution of an electron wave 
packet, interacting with two local spins (located at sites 800 
and 801). For the electron spin we use a tight-binding model 
with nearest neighbour hopping only; for reasons discussed in 
the text we use k — -k 12. For this figure the coupling with 
local spins is given by Jo = 2.0to, and the coupling between 
local spins is set to zero (Ji = 0). The choice Ji = causes 
the time evolution of the electron spin to closely resemble the 
one with a single local spin previously reported in Ref. |25| . 
Subsequent time slices are displaced vertically for clarity. 
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spin interacts with the stationary spins and scatters. If 
there was only one stationary spin, the scattering would 
lead to 4 possibilities for the electron wave packet f25| : it 
can either be reflected or transmitted, with either spin-up 
or spin-down. With two (or more) interacting stationary 
spins, inelastic scattering is also possible. The choice of 
parameters in Fig. [2] is such that the result is similar 
to that expected from a single spin (Ji = here); af- 
ter interacting with the local spins the wave packet both 
reflects and transmits with both spin components. The 
scattering is elastic which means the associated wave vec- 
tors are ±7r/2, so that no spreading of the wave packet 
occurs as time progresses (there is some intrinsic spread 
because two neighbouring scattering sites are involved). 

The 'flnal state' of both the electron and the local spins 
is readily defined by waiting for a period of time after 
which the various electron components have separated a 
reasonable distance from the local spins. This is clear 
from the figure (the latest times shown clearly fulfil the 
above requirement) but we will encounter special param- 
eter regimes where this definition is not so clear, to be 
discussed later. 



(b) 'N' interacting stationary spins 

At the outset we wanted to understand how a (macro- 
scopically) long spin chain interacts with an incoming 
electron spin to understand the effect of a spin current on 
a magnetic layer. With the technology discussed in the 
previous section for treating the time evolution of a cou- 



FIG. 3: (Color online) (a) The ^-component of the electron 
spin long after the electron wave packet has interacted with 
the local spins, as a function of both electron-spin coupling Jo 
and spin-spin interaction Ji for 20 local spins. The outcome 
is sufficiently complicated that we will focus on the problem 
with only two local interacting spins hereafter, (b) Slices are 
plotted as a function of Ji for various values of Jq. As shown 
considerable complexity exists even in these plots. 



pled electron spin/local spins system, the study of rea- 
sonably long spin chains is indeed possible. However, the 
impact on the spin chain is sufficiently complex that this 
program was deemed overly ambitious for the present, 
even if we simply examine the impact on the electron 
spin as it emerges from the spin chain. Looking at 'long 
times' after the interaction, the complexity in a series of 
figures like that in Fig. [2]for various values of Jq and J\ is 
enormous. The summary of such a plot is shown in Fig. 
[21 where the value of the z-component of the electron 
spin is shown after interaction with a spin chain consist- 
ing of 20 coupled S' = 1/2 spins. As a function of the 
interaction parameters Jq and J\ there are quite a num- 
ber of visible ripple-like structures which no doubt are 
related to the excitations that are populated through the 
inelastic scattering channels. This interpretation is rein- 
forced by the observation that, for smaller spin chains, 
the number of ripples is reduced, as the number of possi- 
ble internal excitations is reduced. Slices for fixed values 
of Jo are illustrated in Fig. [Sb, and again it is difficult 
to interpret all the various ripples. For this reason we 
focus, in the rest of this paper, on the simpler system 
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FIG. 5: (Color online) The ^-component of the electron spin 
long after it has interacted with the local spin system, as a 
function of the Kondo coupling Jq. The solid (red) curve 
is the result for a single local spin with S = 1/2 [2^. Note 
that the maximum spin-flip occurs at an intermediate value of 
Jo ~ 2.3to [1^; when two local spins are present the result is 
similar, whether they are non-interacting ( Ji = 0) or strongly 
interacting (Ji = lOto). As one would expect the degree 
to which the incoming electron can reverse its spin is much 
higher when interacting with more than one local spin. 



FIG. 4: (Color online) (a) As in Fig. U the z-component of the 
electron spin long after the wave packet has interacted with 
the local spins, as a function of both electron-spin coupling 
Jo and spin-spin interaction Ji for 2 local spins. This plot is 
discussed extensively in the text. Note the horizontal band 
of strong spin-flip (dark) centred around Jo — 2to, broken 
only near Ji « l.Oto- Smaller Ji values result in independent 
behaviour by the 2 localized spins, while larger values of Ji 
result in strongly coupled behaviour by the 2 local spins. A 
prominent but very slight change occurs along the vertical 
line at Ji = Ito, and a very obvious trough (i.e. a valley as 
far as the z-component of the electron spin is concerned) of 
spin-flip occurs as shown (in dark color) sloping up towards 
the right and exiting the graph at (Ji, Jo) ~ (4fo, Wto)- (b) 
Slices are plotted as a function of Ji for various values of Jo. 
For Jo — 5to,8to there is a definite valley corresponding to 
the dark 'trough' just mentioned in the first plot, while, for 
Jo = 2to, the behaviour is more complicated. 



where there are only two coupled local spins. 



(c) Two interacting stationary spins: inelastic 
scattering 

We first examine the long time behaviour of the elec- 
tron spin. Fig. |4] illustrates (in a color plot) the z- 
component of the electron spin once it has essentially 
left the vicinity of the two local spins, as a function of 
the Kondo coupling between electron spin and local spin, 
Jo, and the coupling between local spins, Ji. Curves are 
shown for the same quantity in Fig. Hl^b), for specific 



values of Jq, as shown; these correspond to horizontal 
sweeps across the first plot. In Fig. [5] vertical sweeps 
across the first plot in Fig. [3] are shown, along with the 
result for a single local spin The sweeps are plot- 

ted for extreme values of Ji, and avoid the complicated 
region characterized by a 'trough' (colored dark) of sig- 
nificant spin flip rising upwards to the right, and leaving 
the plot area at (Ji, Jo) ~ (4, 10)io. This 'trough' region 
will be discussed in detail in the next subsection. 

A considerable amount of information is contained in 
Fig. 21 The horizontal band of strong spin-flip (dark) 
centred around Jo = 2to is further illustrated for specific 
values of Ji in Fig. [5l as a function of Jo (the dark 
horizontal band in Fig. [4^ corresponds to the minima 
visible in Fig. [5|). Whether or not the local spins are 
strongly coupled, the net effect on the electron spin is 
similar, and in qualitative agreement with what happens 
when only a single localized spin is present [1^ (solid 
(red) curve in Fig . [5 | ) As already described for a single 
local spin 
near Jq 



271 130| , the maximum spin flip occurs 
2io; for very small values or very large values 
of Jq the impact on the electron spin goes to zero. 

The reaction of the local spins does depend on the value 
of the coupling between local spins, as illustrated in Fig. 
[SI where the z-component of the two local spins arc shown 
as a function of time for various values of Ji . For zero 
coupling they react independently (except the second lo- 
cal spin 'sees' only part of the incoming electron spin, 
because it has already scattered and spin-flipped off the 
first), while for low coupling some precession occurs. At 
high values of the coupling, the two local spins are essen- 
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FIG. 6: (Color online) The z-component of the two local spins 
for different spin-spin interaction strengths Ji, all for Jo — 
2to. For Ji — the two spins are essentially independent of 
one another, while for Ji = lOto the two local spins are locked 
together with the same value as a function of time. 



tially locked together. 

Referring again to Fig. 31 a subtle change occurs as Ji 
passes through for all values of Jq > ig; this is more 
clearly seen in Fig. [Hb), where a small rise occurs in the 
z-component of the electron spin as Ji/to crosses unity. 
For Jo = 2to the increase is considerable, followed by a 
peak and then a monotonically decaying result. This is 
in contrast to the other two curves which also show a 
minimum. In fact these two curves are more 'generic'; 
inspection of Fig. IH^a) shows that Jq = 2to passes right 
through the middle of the dark band which was discussed 
above. This region of the Jo — Ji phase diagram is fairly 
complicated — the three energy scales are all similar in 
size and no simple picture emerges. 

Focusing on the larger values of Jo, the small increase 
in the z-component of the electron spin shown in Fig. 
IDJb) (also visible in Fig. IDJa) as a faint but abrupt break 
along the vertical line Ji = to) can be understood as 
follows. First note that this increase signals a decrease 
in the spin-flip interaction. Recall that the electron spin 
is propagated with wave vector k — 7r/2. This means 
that its kinetic energy is effectively 2to — the dispersion 
relation e{k) = — 2tocos(fca) gives e{k = tt/2) = 0, but 
2to is the energy with respect to the bottom of the band. 
Thus, the electron has a maximum energy 2tQ that can be 
deposited into the local spin system through the Kondo- 
like coupling Jo. On the other hand, for a two spin system 
there is only one non-zero excitation energy — it is E^x ~ 
2 Ji — and this is essentially the spin wave energy for a 
two spin system, as can be readily ascertained from the 
solution to the problem of two ferromagnetically coupled 
Heisenberg spins 34] . For Ji > to this mode of inelastic 
scattering is no longer possible, so the amount of spin-flip 
scattering decreases, as indicated in the figures. 

An explicit demonstration of this mode of scattering 



FIG. 7: (Color online) A series of snapshots of the electron 
wave packet, with both spin up (solid, red curves) and spin 
down (dashed, green curves). Note that spin down compo- 
nents are scattered both elastically (same speed as incoming 
wave packet) and inelastically (slower speed, indicated by a 
more vertical profile on this plot). The scattering occurs off of 
two local spins, located at sites 800 and 801, ferromagnetically 
coupled with Ji/to = 0.8; we used Jo /to — 2.0. 



is provided in Fig. [3 where a series of snapshots of the 
electron wave packet is shown as a function of position. 
In contrast to Fig. [5] a second set of peaks is evident, 
all in the spin-flip channel (i.e. z component of elec- 
tron spin is —1/2) moving more slowly (hence inelastic 
scattering) both to the left (reflection) and to the right 
(transmission). As Ji to the speed of this wave packet 
approaches zero (so the extra wave packets will appear 
almost vertically in a plot like Fig. [7|). For more and 
more coupled local spins many more inelastic channels 
are available for scattering, which in part explains the 
complexity in Fig. [3l 

(d) Two interacting stationary spins: the 
non-equilibrium bound state (NEBS) 

The most striking feature in Fig. |4]is the trough (dark 
color) that extends upwards to the right, and exits the 
graph at (Ji, Jo) ~ (4, 10)fo- This trough represents a 
domain in the coupling space in which the spin-flip inter- 
action persists more than expected, and is roughly associ- 
ated with a 'resonance' behaviour. The evidence for this 
is very difficult to glean from the numerical calculations 
— we will have more to say based on analytical work to be 
presented in the next section. Nonetheless, examination 
of the numerical results for a particular set of parame- 
ters on a logarithmic scale shows an unusual feature, as 
illustrated in Fig. [8l for relatively high parameter val- 
ues of electron-spin couphng, (Ji, Jo) = (3.1,8)io- On 
this scale the Gaussian wave packets are outside the dis- 
played region at the latest times shown (note that time 
progresses as one moves down from curve to curve, op- 
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FIG. 8: (Color online) A series of snapshots of the electron 
wave packet, with both spin up (solid, red curves) and spin 
down (dashed, green curves). Note that time progresses for- 
ward as one moves from curve to curve downwards, and also 
note the logarithmic scale for the ordinate. By the last times 
shown the usual Gaussian wave packet peaks have disap- 
peared off to the sides; what remains, however, is a small 
peak located near the local spins. We refer to this as a non- 
equilibrium bound state (NEBS); justification for this name 
will come in the next section. Note that this small peak ex- 
ists only in the flipped spin channel. The scattering occurs off 
of two local spins, located at sites 800 and 801, ferromagneti- 
cally coupled with Ji/to = 3.1 and with a Kondo-like coupling 
Jo /to ~ 8.0; with reference to Fig. |4] these parameters place 
us in the middle of the dark colored trough of enhanced spin 
flip scattering. 
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FIG. 9: (Color online) The Sx, Sy, and Sz components of the 
two local spins for the parameter set discussed in the previous 
figure. Note that Sx and Sy remain equal to zero (due to 
the initial conditions, as explained in Ref. H^), while the Sz 
components change, although in reverse order than one would 
naively expect. This is an instance where the classical notion 
of a 'spin vector' that rotates into the direction of the spin 
current while maintaining a constant magnitude is completely 
inapplicable. 



a function of Jq and Ji. We now turn to an analytical 
approach to gain some further insight into the problem. 

ANALYTICAL PLANE WAVE APPROXIMATION 



posite to the progression shown in previous plots). The 
feature in question is the rather small peak located at 
the local spin sites (near site 800 and 801) that persists, 
albeit with strongly diminishing amplitude, for all times 
shown. This peak forms only for the spin down compo- 
nent of the electron; its amplitude decays away in both 
spin channels presumably through a diffusive process, so 
eventually the electron has scattered entirely. We refer 
to this state as a non-equilibrium bound state (NEBS); 
this name will be further justified in the next section. 

In Fig. [5] we show the various components of the local 
spins as a function of time, along with the electron spin. 
The Sx and Sy components remain fixed at zero (because 
of the initial conditions on these spins [1^), while the Sz 
components flip partially and remain at the same value 
long after the flipping process has terminated. In the in- 
termediate stages, however, they are not locked together, 
and remarkably, the 2nd spin flips before the first. This 
reversal of the expected order of flipping occurs only for 
parameters in the 'trough' region; otherwise the local spin 
first encountered by the incoming electron spin is the first 
to flip. While this phenomenon is clearly connected to 
the NEBS, we do not have a simple explanation for the 
spin flip reversal. 

These results illustrate the variety of different be- 
haviour possible for the spin flip scattering process as 



(a) a change of basis 

The problem of an incident spin represented as a plane 
wave scattering off of an impurity with a contact Kondo- 
like spin-spin interaction was solved analytically in Ref. 
[2^. In that problem we made use of the initial condi- 
tions and conservation of angular momentum to simplify 
the problem. Here we do the same, and utilize initial 
conditions such that the Sz component of the incoming 
electron spin is +1/2, while those of the two stationary 
spins are each —1/2. 

The one-dimensional version, written in free space, has 
a Hamiltonian which can be written: 

H = - — — -2jQ[&-Si6{x)+a-S2S{x~a)]~2JiSi-S2. 
Zm ax'^ 

(11) 

The wave function for this problem consists of a spatial 
component which describes the electron spin amplitude, 
and a spin part which describes the spin state of the 
incoming electron and the two stationary spins (here lo- 
cated at positions x — Q and x = a). The Hilbert space 
concerning the spin degrees of freedom has an overall size 
of 2'^ = 8 (for 5=1/2 spins) However, utilizing the con- 
servation of total Sz reduces this number to 3. As already 
stated, the initial state, in Dirac notation, is | tii); where 
the first arrow represents the 2;-component of the electron 
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spin, and the next two arrows indicate the respective z- 
components of the two local spins. Once the electron 
spin interacts with the local spins, two more spin states 
are possible, | |ti) and | lit)- In our numerical results, 
we followed two separate routes: in cases with the initial 
configuration as depicted here, we used this fact to reduce 
the Hilbert space to these three spin states, which sped 
up the calculations considerably. Alternatively, when the 
initial configuration was not so straightforward (and did 
not have a definite total S'z, for example), we used all 8 
basis states. 

When we begin with an initial configuration like | tii), 
we can combine these spin states into combinations with 
both good total Sz and good total S to give rise to the 
following basis set [sst : 



f 



l^i) = ^(1 iiT> + UTi> + l Tii» 



(12) 



m = ^(iuT> + im>-2iTU)) (13) 



IV'3) = -^(iuT)-im»- 



Writing the wave function as 



(14) 



|V(x)) = /i(x)IV'i) + j[x)\^2) + g{x)\i,^), (15) 



then appropriate projection on to the spin basis states 
results in the three equations: 



-2Jo-(5(x)+^(.x-a))/. 
+ ^oy ['5(x)(2/ - V3g) + b{x - a) (2/ + V3.g)] 
-|^0-V3J„y[5(.)-^(.-a)]/ 



<^el.f 

{eel - 2Ji?i^)g 



(16) 
(17) 
(18) 



r 



where Cei is the kinetic energy of the incoming electron. 
Note that the first equation results from the S'tot = 3?i/2 
sector, and remains decoupled, while the second two are 
part of the 5tot = 5'totz = — ?i/2 doublet. The state 
with spatial wave function g{x), governed primarily by 
the third equation, exists exclusively because of the pos- 
sible inelastic scattering process. Eq. ([T^ indicates that 
it contains only the spin down component of the scat- 
tered electron, and, given our initial conditions, exists 
only after scattering. It is 'fueled' through the f{x) com- 
ponent, which, as Eq. ()13|1 indicates, contains a compo- 
nent corresponding to the incoming electron spin (with 
az = +h/2). That the g{x) component corresponds to 
inelastic scattering is indicated by the eigenvalue on the 
right-hand-side of Eq. (fT8|l . with value eei — 2h^Ji, which 
shows that an energy 2h^Ji is left behind in the form of 
a spin wave excitation in the local spin system, as ex- 
plained in the previous section. The first two equations, 
Eqs. (fT6|) and p?)) . each have eigenvalue ed, showing 
that the kinetic energy of the incoming electron is con- 
served (elastic scattering). Note that this still results in 
spin-flip scattering; it is just that the two local spins are 
scattered by the same amount, so that they remain in 
their coupled ground state. 



(b) a re-examination of the numerical solutions 

Eqs. (jlGHlSp can be readily solved analytically, and 
we will come to that solution shortly. However, already 
Eqs. (|12m4p serve the important task of directing our at- 
tention to specific linear combinations of the spin states, 
as indicated. The numerical solutions presented in the 
previous section were classified only according to the z- 
component of the electron spin. We now essentially re- 
plot those results, as separate amplitudes h{x), f{x), and 
g{x), in Fig. [TUl Note that Eqs. (|16m8p were derived 
for the continuum model defined by Eq. (jlip : nonethe- 
less the role of the various amplitudes, described at the 
end of the previous subsection, applies equally well to the 
numerical results of the original tight-binding model. 

To demonstrate this, in Fig. [TOTa') we plot the magni- 
tude |ft,(a;)p vs. position for a number of time slices, for 
three different values of Ji. As predicted by Eq. (flBjl. 
there is no dependence on Ji. It is important to note 
that these results still represent numerical solutions to 
the tight-binding model presented in the previous sec- 
tion; while we could have used the spin components as 
listed in Eqs. (|13lll4p as a basis set, these numerical solu- 
tions do not 'use' the analytical structure of Eqs. (I16m8|) . 
Hence only one set of curves is visible (for Ji = lAto) as 
this set is identical to and covers entirely the sets corre- 
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FIG. 10: (Color online) The time evolution of the magnitudes 
(a) |/i(a::)p, (b) \ f{x)\'^, and (c) [(/(s)^, as defined by the basis 
set, Eqs. (|12I14|I . These plots apply for Jo = 2to and the 
values of Ji indicated. Note that in (a) the plots are identical 
for all 3 values of Ji, as motivated by the structure of Eq. 
(|16l) . In (b) and (c) differences are apparent; note that in (c) 
no amplitude is present before the time of scattering, and, 
furthermore, as one enters the 'trough' region (Ji = lAto) 
\g{x)f consists of a single sharp peak near the local spins. In 
time this peak diffuses outwards, but there is no wave packet 
component. 



spending to the other two values of Ji . 

In contrast, the other two components, plotted in Figs. 
[TUl (b) and (c), are dependent on the value of Ji. In 
both cases the amplitudes of transmitted and reflected 
wave packet depend quantitatively on the value of Ji. 
Note, moreover, that the amplitude g{x) has no 'incom- 
ing' wave packet. As explained earlier this amplitude is 
generated entirely by the scattering process. Also note 
that for Ji/to < 1 (i.e. Ji/to = 0.8 in Fig. [TOl) the 
slow moving piece belongs entirely to g{x) while the fast 
moving one belongs entirely to f{x). 
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FIG. 11: (Color online) A plot of \g{x)f (see Eq. ^) vs. 
X. This 'snapshot' is taken immediately following the initial 
scattering of the electron spin with the two local spins, situ- 
ated at sites 790 and 810, i.e. they have been separated for 
clarity. We use Ji = S.lto, so it is clear that for parame- 
ter values that fall on the 'trough' (Jo = 8to) the component 
of the wave function associated with inelastic scattering (i.e. 
g{x)) is significantly enhanced (almost two orders of magni- 
tude) compared with the region away from the trough. 



To see the role of the g component of the state more 
clearly, we separate the two local spins by 20 sites, and 
project out the g component from the numerical solution, 
using Eq. (fT5|) . In Fig. [Tl] we show on a log scale the 
magnitude of the g component, |.g(a:)p as a function of 
position; the two local spins are now located at sites 790 
and 810. The parameters {Jo/to, Ji/to) — (8,3.1) (solid 
curve) situate the regime on the 'trough' so apparent in 
Fig. m whereas {Jo/tQ, Ji/to) = (2,3.1) (dashed curve) 
puts one well away from the trough. This snapshot is 
taken immediately after the initial scattering takes place, 
and it is clear that the g component is almost 2 orders 
of magnitude larger on the trough (solid curve) than off 
(dashed curve). A similar result holds for large values of 
Jo. 



(c) analytical solution 

An analytical solution of the problem with plane waves 
through Eqs. (|16m8p is possible, though tedious. One de- 
fines three regions in space, and defines the wave function 
in a piecewise continuous manner, as is done commonly 
in undergraduate physics texts. 

With k = 



' 2mec\/h , and Q — 
the wave function can be written 



h{x) 



fine'''- 



rue 



-ikx 



X <Q 
<x < a 
X > a 



X <Q 
< X < a 
X > a 



eei)/r, 



(19) 



(20) 
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and 

gix) ^ < giie-Q"" + siie'^'=, 0<x<a \. (21) 
[ giiie^'^'', x> a ] 

Four conditions relate the various coefficients defining 
h{x) to the incoming amplitude hj in Eq. (fT9|) : simi- 
larly eight conditions determine the / and g coefficients 
in terms of the incoming amplitude //. The results for 
h{x) are standard and can be found in many undergrad- 
uate texts, while, for / and g, the result is not stan- 
dard, but is nonetheless straightforward. Also note that 
we have written the wave function for the more relevant 
condition 2h? Ji > ed, in which case the function g{x) 
is exponentially decaying; the alternative 2fi^Ji < ed is 
straightforward and gives a propagating wave solution, 

with wave vector q = ^J2m{tc\ — 2h^Ji)/h^. This latter 
case corresponds to the situation whereby a spin wave 
excitation is energetically allowed, so that a spin-flipped 
wave packet will emerge from the stationary spins at a 
reduced speed, as we have already seen in the numerical 
solution in Fig. [T] 

When inelastic scattering is not allowed by energy con- 
siderations, it is not clear what will happen. Our intu- 
ition suggests that the stationary spins will respond as 
one, and so the spin-flip process will resemble that ex- 
pected for scattering from a single spin (which, as we 
demonstrated earlier, is not so different from scattering 
off decoupled stationary spins (Ji = 0)). Inspection of 
Fig. d] shows that this is indeed the case, except for the 
trough region previously identified. It is precisely in this 
regime that a peculiar enhancement of spin-flip scattering 
occurs, which we now argue is connected to the effective 
bound state (NEBS) defined by Eq. (|2T|) . 

The solutions can readily be written down by using 
the definitions, a = Jo/k and (3 = Jq/Q (the mass m 

is set equal to unity), and the terms v = a^^l — ^(1 — 

g-Qagjfca)| and u = a|l - ^{1 - g-^'^e-''''"-)^ appear 

often. Note that for real Q these are complex conjugates 
of one another. However, these expressions (and the ones 
immediately following) are valid for high electron kinetic 
energy as well, where ed > 2Ji, and so it follows that 
Q = —iq with q now real, and u and v are no longer 
complex conjugates of one another. We find, for example, 

gu ^ VI J_±iv-iu^^ 
fi 2 ' {1 + ivy + u2e2*fe'^ ' ^ 

with similar expressions for the other coefficients. To see 
how effective the spin flip process is, we can calculate ei- 
ther the expectation value of the electron spin, (uz), or 
the spin torque, iVzx [IS]- For the two local spin system 
used here, the latter is given in terms of the former as 
-^zx = k{l/2 — ((Tz)). As in the earlier numerical results. 




FIG. 12; (Color online) (a) Plot of the expectation value of 
the z-component of the electron spin, (ctz), as a function of the 
two coupling parameters, Ji and Jo, based on the plane wave 
solutions to this problem. The range of Ji and Jo is consider- 
ably extended compared to Fig. |4]to emphasize the presence 
of the trough (shown in dark color) that extends upwards and 
to the right. In (b) we show a plot for 1/2 — 0.2 * |<7//P/1//P 
for the same parameters; the trough is immediately identi- 
fiable in this plot, which reinforces our contention that this 
region of enhanced spin flip scattering is associated with the 
non-equilibrium bound state (NEBS) represented by gu (a 
plot of Si I yields similar results) . 



the quantity (az) will remain near 0.5 (the initial elec- 
tron spin value) if very little spin flip occurs, whereas 
this quantity will deviate most from 0.5 (and even be- 
come negative) when significant spin flip occurs. Note 
that with the plane wave solution given in Eqs. (|19ll2ip . 
the problem is no longer time dependent; one envisions 
a continual influx of current (this is //) while reflected 
and transmitted plane waves (of both spin type) take on 
'steady-state' values 

The calculation of (ctz) is straightforward; we use an 
integration region —L < x < +L, and we allow L oo. 
The plane wave regions outside the local spin region then 
dominate, and, for real values of Q, we obtain 

(^.) = Y^{5 - ^V2Re{hinfhi + ^ir})}, (23) 

while, for pure imaginary values of Q, the expression for 
(ctz) is somewhat more complicated. 

In Fig. [T2]we show (a) (ctz) and (b) \gii\'^ as a function 
of Ji and Jo to emphasize the connection between the re- 
gion (described as a 'trough') of enhanced spin flip scat- 
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tering and the non-equihbrium bound state (NEBS). The 
range of both Ji and Jq is considerably extended com- 
pared with Fig. m nonetheless the qualitative similarities 
are striking; clearly the analytical solution captures the 
essence of the numerical one. Furthermore, the analytical 
approach has allowed us to make the association of the 
trough of enhanced spin flip scattering with the NEBS. 
Quantitative details differ, in part because the numerical 
results are based on a tight-binding model whereas the 
analytical ones utilize a quadratic dispersion for the itin- 
erant spin. A specific example is given in Fig. 1131 where 
both {(Tz) and \gii\^ are plotted as a function of Jq (for 
a specific value of k and Ji). Clearly the peaked region 
in \gii\^ (near Jq ~ 15) corresponds to the dip in (ctz), 
showing strong evidence for the role of the NEBS in en- 
hanced spin flip scattering. For large values of Ji >> ed, 
Eq. ((22l) simplifies somewhat; we get 



■ 2 3 / Jo\ sin fca -I- (cos fca -I- 2u sin fca) , , 
^A\q) l + 4i;2(cosfca-f wsinfca)2 ' ^ ' 

where v ~ ^ [l ~ I^]- Similar analytical expressions 
can be readily attained for all the coefficients, but they 
are of limited value. 

Eq. ([24|) is also plotted in Fig. [131 where it is seen 
to be very accurate (in fact, it is fairly accurate all the 
way down to Ji « 2). The peak region in \gii\'^ follows 
roughly a dispersion relation 



Ji 



2 



64 



' 



(25) 



and, as has been emphasized already, this corresponds to 
the region of most intense spin flip scattering (the 'dark 
trough' region of previous figures). Thus, when Jq and 
Ji are tuned to satisfy Eq. ([25]) we find an enhanced spin 
flip process. 




FIG. 13: (Color online) (a) Plot of < > (solid (red) curve) 
vs. Jo for specific values of Ji and k, as indicated. Also shown 
is the coefficient \gii'^ (solid (green) curve), which shows a 
peak precisely where < ct^ > has a significant dip, indicative 
of enhanced spin flip scattering. Also shown (symbols) is the 
result of an approximate analytical expression derived in the 
text. Agreement is extremely good. 



and Ji. As is clearly evident in (a) and (b), the trans- 
mission of both spin species is noticeably enhanced in 
the trough region. Coincidentally the spin-up reflected 
component is decreased, while the spin-down reflected 
component shows an increase. The increase in the trans- 
mitted spin-up component of the electron is not through 
'direct' transmission. Rather it is achieved through the 
spin-flip interaction that generates the component with 
amplitude g discussed earlier in this section. Recall that 
in this parameter regime this g-component does not exist 
outside the local spins; it first transforms into the compo- 
nent with amplitude /, which represents a propagating 
wave with both spin- up and spin-down species. These 
plots therefore reinforce the idea that the electron goes 
through a two step 'virtual' spin flip interaction (creation 
of the NEBS) in the trough region. 



(d) transmission and reflection amplitudes from the 
numerical solutions 

Having established the idea of a NEBS we once again 
re-examine the numerical solutions. In particular, one 
important property from the experimental point of view 
is that the stationary spins can act as a spin barrier. We 
have already shown that a large electron-spin interaction 
(Jo) works as a high potential barrier for the incoming 
spin. In our framework, for instance, a large electron- 
spin interaction acts to prevent any spin-up component 
of the electron from transmitting through the stationary 
spin system. However, in the Jq — Ji phase space, at the 
onset of the trough described, for example, in Fig. [4l the 
transmitted component of the spin-up (and spin-down) 
electron is enhanced considerably; this is illustrated in 
the 4 plots shown in Fig. [TH where both transmitted 
and reflected intensities are plotted as a function of Jq 



CONCLUSIONS 

We have modeled spin current-induced spin torque in 
the quantum regime with a lattice, on which an itinerant 
spin (constructed as a wave packet) moves with a kinetic 
energy given by a tight-binding dispersion, to represent 
the spin current. Any number of ferromagnetically cou- 
pled spins can then be flipped by repeating the process 
described here with more itinerant electrons, i.e. a cur- 
rent. As described in Ref. [H, this then requires a density 
matrix description. We have focussed on just two coupled 
local spins, since this small system contains the essence of 
the processes we believe are responsible for spin torque: 
(i) direct spin flip without internal excitation of the local 
spin system, and (ii) spin flip through inelastic scatter- 
ing, either real or virtual. The first process exists even 
for a single local spin, and has been explored previously 
by us. The second process is the primary subject of this 
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FIG. 14: (Color online) The (a) transmitted spin-up, (b) 
transmitted spin-down, (c) reflected spin-up, and (d) reflected 
spin-down intensities as a function of Jo and Ji . These results 
are obtained with the same conditions as in Fig. |4l See the 
text for further description. 



paper, particularly in the regime where, energetically, the 
itinerant spin becomes momentarily bound in the local 
system, a phenomenon which we have called the non- 
equilibrium bound state (NEBS). The description here is 
for a one dimensional system, but the NEBS should also 
be present in three dimensions. 

An analytical plane wave approach, using a parabolic 
dispersion for the itinerant spin, helps to elucidate the 
nature of the spin flip processes. A scattering channel 
through which a local spin singlet is generated is respon- 
sible for the enhanced spin flip scattering along a 'trough' 
in the (Jq, Ji) phase diagram. This trough is reasonably 
well described in the plane wave approach by the relation 
Ji = 9 J2/64 + eel/2. 

An experimental observation of the NEBS would be 
straightforward provided at least one of the parameters 
Jo, Ji, or Eel can be tuned in a particular system. In 
this way the probability of spin flip can be monitored as 
a function of parameter space, and the NEBS would be 
identified by a well defined region of enhanced spin flip, 
corresponding to the 'trough' in Fig. [H 

One interesting consequence of our calculation is the 
possibility of using the spin chain as an effective spin 
filter. By tuning the parameters to correspond to the 
regimes of enhanced spin flipping, the spin up electrons 
will be flipped while the spin down ones will be unaf- 
fected. This effect can be achieved not only for the 2-spin 
chain but also for the longer chains, as shown in Fig. [S] 
The resonant 'trough' provides a controllable spin filter 
through the inter-spin coupling Ji . 
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